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ABSTRACT 


An  overall  review  of  studies  of  electrostatic  interactions  with 
fluid  flows  is  presented.  Review  of  prior  studies  involving  ions  with 
free  convection,  and  channel  flows  is  presented  along  with  prior  studies 
of  field  effects  on  condensation.  The  cause  of  large  changes  in  laminar 
channel  flow  with  ionization  was  studied  both  analytically  and  experi¬ 
mentally.  The  changes  were  determined  to  result  fran  an  electrically 
induced  destabilization  of  the  flow.  External  ion-flow  boundary  layer 
actions  were  studied  and  sane  influence  on  transition  was  found.  Sig¬ 
nificant  changes  to  frost  formation  and  flame  heat  transfer  under  the 
action  of  an  impressed  E  field  were  found.  The  possible  use  of  ions 
injected  into  a  gas  stream  as  a  flow  diagnostic  technique  was  explored 
and  is  considered  a  very  logical  area  for  future  exploration. 
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I.  INTRODUCTION 


This  report  presents  a  broad  summary  of  the  work  that  was  con¬ 
ducted  under  the  project  sponsored  by  the  U.  S.  Army  Rest  ch  Office  - 
Durham.  The  specific  area  of  investigation  is  the  interaction  of 
electrostatic  fields  on  fluid  flows  and  boundary  layers.  8ince  a 
broad  range  of  interactions  is  possible ,  the  scope  of  this  report  is 
wide,  ranging  from  laminar  channel  flows  through  ion-diagnostic  tech¬ 
niques  . 

The  report  itself  presents  same  background  material  accomplished 
by  the  principal  investigator  prior  to  this  program  to  indicate  the 
overall  continuity  of  effort  in  this  particular  field  of  endeavor. 

The  report  consists  of  several  sections.  Prior  work  is  discussed  is 
the  sections  on  free  convection,  dielectrophoresis,  condensation,  and 
the  initial  work  on  channel  flows.  The  more  complete  reports  or  papers 
summarizing  each  of  these  areas  are  listed  in  the  references.  The 
remaining  work  discussed  herein  was  accomplished  under  this  project 
with  the  exception  of  a  portion  of  the  ion-flow  diagnostic  tests.  The 
use  of  ions  as  a  diagnostic  tool  in  the  channel  flow  tests  was  done 
under  this  program,  as  well  as  the  study  and  plain..  <13  of  the  use  of 
ion-flow  diagnostics  in  external  flew.  The  actual  preliminary  testa 
of  the  external  flow  diagnostics  was  accomplished  by  the  principal 
investigator  at  the  U.  S.  Army  Aeronautical  Research  Laboratory, 

Moffett  Field,  California. 

Specific  work  conducted  under  this  program  included  an  evaluation 
of  the  ion-channel  flow  interaction  found,  study  of  electrostatically 
induced  secondary  flows  in  channels,  ion  effects  on  external  boundary 
layers  using  an  interfereameter  technique  and  also  a  traversing  probe 
technique,  electrostatic  effects  on  flip-flop  fluidic  devices,  effects 
of  fields  on  flame  heat  transfer,  and  effects  of  fields  on  frost  forma¬ 
tion.  The  complete  reports  prepared  on  the  above  topics  are  indicated 
in  the  list  of  references. 


II.  BACKGROUND 


A.  TYPES  OF  ELECTROSTATIC  FORCES 

The  study  of  the  effects  of  electrostatic  fields  on  fluid  flcsrs 
centers  on  the  possible  ways  in  which  electrical-fluid  interactions 
can  exist.  Two  broad  categories  can  be  delineated,  charge  effects  and 
nonuniform  or  gradient  influences.  In  the  first  category  fall  all  the 
usual  coulomb  force  effects  which  depend  on  the  existence  of  a  free 
charge,  and  in  the  second  category  are  more  subtle  influences  in  which 
gradients  of  properties  of%  the  fluid  medium  or  the  field  are  important. 
If  one  considers  the  body  forces  which  can  arise  due  to  electrostatic 
effects,  then  the  following  equation  depicts  the  body  force*: 
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Pc  ®ext  +  \  6*®* 


(1) 


The  charge  force  is  the  first  term  on  the  right  hand  side,  and  the  non- 
uniform  effects  are  given  by  the  remaining  two  terms.  In  the  case  of 
a  gas,  the  above  relation  assumes  the  more  familar  form: 

fgu  *  Pc  e"4  **  <*> 


The  interesting  term  in  Eq.  (2)  is  the  grad  E2.  A  similar  term  is 
found  as  an  approximation  in  the  case  of  liquids.  The  grad  I2  terns 
implies  that  a  body  force  can  exist  on  a  electrically  neutral  fluid 
particle  if  a  nonuniform  electrical  field  exists.  The  other  category 
of  forces  involved  in  electrostatics  are  the  surface  stresses  due  to 
a  field  and  are  of  the  form: 
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(3) 


The  nature  of  the  charge  involved  in  the  term  p^E  of  Eq.  (1) 
(wheie  pc  is  the  charge  density)  depends  upon  the  medium  as  well  as 
the  means  of  generating  the  charge.  That  is,  in  an  ordinary  gas  at 
low  temperatures  (below  1000°F)  the  ions  will  consist  of  those  pro¬ 
duced  by  external  radiation,  or  impressed  fields.  They  will  Include 
Qs  and  N2  ions  in  ordinary  air  as  well  as  electrons.  Existent  impuri¬ 
ties  will  also  become  charged  and  act  as  massive  ions.  In  flames  a 
host  of  nonequilibrium  combustion  ions  are  present  including  such 
specie  as  HsO*.  In  liquids  the  charged  particles  of  Importance  to 
these  studies  generally  are  the  impurities  which  exist  as  contaminants 
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in  the  liquid  medium.  Such  impurities,  (minute  particles)  can  acquire 
a  charge  either  because  of  induced  ionization  by  an  external  field  or 
through  the  natural  double  layer  charging  process  familar  to  electro¬ 
chemistry.  When  considering  the  charge  effects  on  fluids  one  of  the 
singularly  most  important  factors  is  the  mobility  of  the  ions  involved. 
Mobility  is  a  measure  of  the  rate  at  which  a  charged  particle  will  move 
through  a  fluid  under  the  action  of  an  impressed,  field.  Typically  an 
electron  has  an  extremely  high  mobility,  whereas  a  charged  colloid 
particle  has  a  very  low  mobility.  The  mobility  also  serves  as  an  indi¬ 
cator  of  the  ability  of  a  charged  particle  to  couple  with  the  fluid 
medium.  An  electron  does  not  couple  well  with  the  fluid  whereas  the 
colloid  ion  does.  Hence  one  can  expect  the  interactions  of  electric 
fields  with  the  fluids  to  be  strongly  dependent  on  the  mobilities  of 
the  particles  involved. 
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III.  POSSIBLE  INTERACTIONS 


Based  upon  the  nature  of  the  electrostatic  forces  which  can  be 
exerted  it  is  possible  to  hypothesize  various  fluid  phenomena  which 
could  be  affected  by  electrostatic  fields  or  charges.  The  actions 
include  the  surfaces  or  interfaces  of  various  liquids,  change  of  phase 
such  as  boiling  and  condensation,  boundary  layer  flows,  convective 
heat  transfer,  flow  attachment,  and  fluid  orientation.  Considerable 
study  has  been  accomplished  in  the  area  of  fluid  surface  effects  by 
Melcher  wherein  he  has  shown  that  both  liquid  surfaces  and  streams  can 
be  destabilized  using  electrostatic  fields.2  The  other  possible  inter- 
actions  will  be  considered  in  greater  detail  in  this  report. 

A  particularly  fascinating  possibility  which  arises  with  these 
interactions  is  their  use  in  fluid  flow  diagnostics.  Since  the  ions 
have  relatively  little  inertia,  it  is  considered  likely  that  they  can 
be  used  as  tracers  to  follow  the  fluid  particles  in  the  stream.  Illus¬ 
trations  of  the  use  of  the  ions  in  this  fashion  also  are  given  later 
in  this  report. 


A.  EXPERIMENTAL  FINDINGS 

The  electrostatic- fluid  interactions  which  have  been  studied  are 
grouped  in  areas  of  free  convection,  internal  flows,  external  flows, 
flow  attachment,  change  of  phase,  flame  heat  transfer,  and  ion- flow 
diagnostics. 


1.  Free  Convection 

The  initial  tests  conducted  to  explore  the  possible  use  of 
the  electrostatic  fields  on  fluid  flows  were  directed  at  studying  the 
influences  on  free  convection  from  a  flat  plate.3  The  test  arrsnge- 
ment  consisted  of  mounting  an  isothermally  heated  flat  plate  in  a 
Mach  Zehnder  interferometer  and  stuyding  the  effects  of  various  applied 
electrical  actions.  The  particular  arrangement  was  chosen  to  use  a 
very  sensitive  instrument  to  observe  small  changes  to  a  very  sensitive 
phenomena,  free  convection.  The  arrangement  is  shown  in  Fig.  1,  where 
the  heated  plate  is  shown  along  with  a  typical  electrode.  In  this 
case  the  electrode  was  a  single  wire  mounted  horizontally  at  the  center 
of  the  plate.  Voltage  was  applied  between  the  wire  and  the  plate.  In 
most  runs  the  electrode  was  positive,  but  runs  were  also  made  with  the 
heated  plate  positive.  Figure  2  illustrates  the  action  of  the  field 
on  the  thermal  free  convection  boundary  layer.  It  can  readily  be  seen 
that  as  the  corona  current  begins  to  flow,  severe  distortion  to  the 
thermal  boundary  layer  occurs.  At  high  currents  the  pattern  is  symmet¬ 
rical  about  the  centerline  directly  under  the  wire.  It  should  be  noted 
in  each  interferogrem  that  the  fine  0.004-incb  corona  wire  is  located 
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Plate  and  wire  electrode  schematic  illustrating  coordinates 


H  4 


at  the  right  end  of  the  support  rods  which  protrude  from  the  lift  side 
of  the  interferogram.  The  wire  direction  is  perpendicular  to  the  plane 
of  the  picture.  On  the  right  hand  side  the  straight  vertical  dark  sur¬ 
face  is  the  heated  plate  and  the  interference  fringes  exist  in  front  of 
the  plate. 

The  severe  distortions  which  were  observed  were  cc..sidcred 
to  be  due  to  the  electric  wind  induced  by  the  corona  discharge.  To 
test  this  hypothesis,  runs  were  made  with  the  plate  in  a  horizontal 
position  with  the  heated  side  facing  down.  Figure  3  illustrates  the 
results.  The  wire  is  horizontal  and  is  at  the  centerline  2  cm  below 
the  plate.  The  ”T"  shaped  rod  provides  reference  distance  for  data 
reduction.  From  this  figure  it  can  be  seen  that  the  severe  distortion 
previously  found  is  absent  and  a  symmetric  well-b^hav  pattern  exists. 
It  was  concluded  that  the  previous  distortion  was  caused  by  the  inter¬ 
action  of  heated  air  up  the  plate  wit  the  corona  wind  flowing  to  the 
plate. 


An  analysis  was  conducted  of  the  drw  ard  :  icing  case,  assum¬ 
ing  the  action  was  due  to  the  electric  w4  d.  Jsinr  he  -lectrostatic 
equations  for  the  electr  .  vind,  and  tv*  moo  n.um  •  trjal  method  for 

the  beat  transfer  analysis,  explicit  solitior  was  obtained. 
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where  Z  *  x/a  and  ra  is  the  wire  radius.  The  results  are  shown  in 
Fig.  4  where  reasonable  agreement  between  theory  and  experiment  is 
shown.  The  primary  conclusion  from  this  work  is  that  the  electrostatic 
actions  are  predictable  and  controllable. 
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HORIZONTAL  PLATE ,  2  Cm  SPACINC 
•ITH  APPLIED  FIELD;  5.0  kv, 

1.05  MICROAMP. 


HORIZONTALLY  MOUNTED  PLATE;  NO 
APPLIED  FIELD,  0.004  INCHES  CORONA 
VIDE  MOUNTED  2  cm  BELOV  CENTER¬ 
LINE  OF  PLATE 


Fig.  3  -  Interferometer  picture  with  horizontal  plate 


Subsequent  work  in  the  area  of  free  convection  effects  has 
been  undertaken  at  the  Air  Force  Institute  of  Technology  and  has  been 
directed  ■‘owards  the  actions  of  electric  wind  at  low  pressures  and 
forced  convective  cells  induced  by  electric  wind.4*® 

2.  Dielectrophoresis 

An  exploratory  study  was  made  of  the  grad  E2  effects  using  a 
simulated  zero  gravity  model.3  Droplets  of  water  were  suspended  at  the 
interface  of  two  nonconducting  liquids.  The  arrangement  is  shown  in 
Fig.  5.  The  central  electrode  just  protrudes  into  the  top  of  the  liquid. 
The  ground  electrode  is  on  the  outside  of  the  container.  After  the 
field  (4,000  volts  a.c.)  is  applied  it  can  be  seen  that  the  final  posi¬ 
tion  of  the  water  droplets  is  at  the  center  electrode  (Fig.  6).  Thus 
the  potential  utility  of  the  nonuniform  field  term  is  demonstrated. 

Work  utilizing  this  principle  has  been  carried  out  by  various  groups  in 
an  attempt  to  provide  for  "zero-G"  fluid  orientation. 

3.  Channel  Flow 

A  fundamental  aspect  of  a  study  of  fluid  flows  is  the  influence 
of  the  viscous  effects,  leading,  for  example,  to  boundary  layer  flows. 
Because  of  the  critical  nature  of  boundary  layers  to  practical  fluid 
mechanic  problems,  study  was  initiated  on  the  possible  effects  of  elec¬ 
trostatic  actions  on  viscosity  and  viscous  flows.  A  simple  experiment 
was  conducted  using  channel  gas  flow  in  which  ions  were  caused  to  flow 
transverse  to  the  main  stream  in  a  radial  direction.®  The  arrangement 
is  shown  in  Fig.  7.  Gas  is  passed  into  the  plenum  at  the  left,  passes 
through  the  round  channel  and  exhausts  at  the  right.  Pressure  taps  are 
located  to  measure  the  pressure  drop  in  the  channel.  The  ions  are 
obtained  by  impressing  a  suitable  voltage  between  the  wire  and  the 
cylinder  wall.  Typical  ion  densities  are  1  x  1010  ions  per  cubic  centi¬ 
meter  in  air  at  N.T.P.  Results  of  the  tests  are  presented  in  Fig.  8 
where  a  typical  friction  factor-Reynolds  number  plot  is  presented. 
Although  study  of  the  ion-molecule  interactions  and  the  Navier  Stokes 
equations  had  indicated  only  a  very  small  influence  of  viscosity,  the 
test  results  show  that  the  friction  factor  has  doubled  in  the  laminar 
flow  regime.  The  velocity  profile  associated  with  this  action  is  shown 
in  Fig.  9.  At  low  Reynolds  numbers  the  profile  changes  from  a  "modified 
parabolic"  shape  due  to  the  presence  of  the  corona  wire  and  no  ion  flow 
to  a  profile  which  is  similar  to  slug  flow.  It  can  be  seen  from  a  study 
of  both  figures  that  the  ion  influence  tends  to  disappear  in  the  turbu¬ 
lent  regime.  Heat  transfer  tests  indicated  an  increase  which  corre¬ 
sponded  to  the  increase  in  frict'on  factor.  The  results  implied  an 
increase  in  effective  viscosity  similar  to  that  reported  for  liquids 
with  applied  fields.7 »a 

4.  Channel  Flow-Approximate  Analysis 

If  direct  solution  of  all  of  the  electrical  field,  current 
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CROSS  SECTION  OF  TUBE 

Fig*  7  “  Experimental  arrangement  of  tube  flow  equipnent 
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and  flow  equations  is  attempted,  it  leads  to  a  complex  nonlinear  equa¬ 
tion.  Another  approach  based  upon  the  actual  experimental  problem  is 
possible.  It  can  be  assumed  that  the  low  velocity  flow  does  not  affect 
the  mechanism  of  the  corona  discharge  materially  and  consequently  the 
corona  charge  density  can  be  treated  independently  from  the  flow.  In 
addition,  because  of  the  length  of  the  channel  used  in  the  tests,  it 
can  also  be  assumed  that  the  transverse  charge  density  distribution  is 
essentially  constant  along  the  axis  of  the  channel.  Figure  10  illus¬ 
trates  the  coordinates  used. 

The  current  equation  along  the  axis  of  the  channel  can  be 
written  as : 


Jz  -  pc  w  ♦  pc  K  Ez 


(4) 


Die  boundary  conditions  on  the  problem  require  that  the  current  path 
must  close  in  the  eystem  and  that  all  charge  that  leaves  the  wire  must 
reach  the  channel  wall.  This  condition  requires  that  no  net  current 
flows  out  of  the  channel,  iz  ■  0  at  the  end  of  the  channel.  Hence 


j[jzdA.O  (5) 

For  the  cylindrical  case  the  radial  current  density  is  given 
by 


Jr  •  ir/2  n  r  1 


(6) 


Using  dlv  J  -  0,  for  cylindrical  coordinates, 


1  dr  Jp  i 

?-sr-  +  "5T 


0 


and  substituting  Jr  from  Eq.  (14) 


dz 


0 


since  ir  must  be  a  constant. 


(7) 


(8) 
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Substituting  Eq.  (4)  into  Eq.  (5)  results  in 


£  (pc  w  ♦  Pc  K  E*)  ^  ■  0 
or 

X  i  E‘  “  ■  -  iL  *  z  “ 


(9) 

(10) 


It  can  be  recognized  that  each  integral  in  Eq.  (10)  represents  the  Bean 
value  of  a  function.  Hence  we  can  write 

/  „  \  (Pcw)m  /..x 

(Pc®z)m  "  “  £  (H) 

and  noting  that  the  electrostatic  body  force  along  the  axis  of  the 
channel  is 


Ft  ■  PcEz 


(12) 


Then  the  mean  value  of  pcEz  is  equal  to  the  mean  value  of  the  body 
force.  Substituting  into  the  fluid  flow  equation  for  channel  flow  with 
a  body  force  one  obtains 


M  vSr  -  j  (pc*)a 


d£ 

dz 


03) 


Non-dimensionlizing  this  equation  for  a  round  tube, 


dSr  A  1  dw  (Pcv)a  w 
— r  ♦ - "  nD 

dTl  dTl  Pcja 


(14) 


where 


Pc*  L2/u  K,  pc* 


Pc 

■  ■■  •  i 

Pcm 


n  .  r/R,  C  -  z/R,  V. 


The  charge  number,  N-,  is  a  new  dimensionless  parameter  which  char¬ 
acterizes  flow  with  ionization.  It  represents  the  ratio  of  the  charge 
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forces  to  the  viacoua  force  in  the  flow.  Assuming  the  charge  density 
ia  constant  one  can  solve  the  equation  to  get 


¥m  ■  ‘ 


(15) 


and  the  friction  factor  becomes 


fD 


(16) 


As  a  typical  case  consider  the  following  values,  K  ■  2  x  IQ**4 
m2/volt-sec,  R  ■  O.OI587  m,  i  ■  10”3  anq>,  channel  length,  i  ■  1.7  meters, 
(i  -  1.783  x  10“®  kg/m-sec,  then  *  2.5  x  10s  charges/c. c.  The  charge 
number  becomes 


■  12. 


This  la  a  surprisingly  large  value.  Substituting  this  value  of  the 
charge  lumber  into  Eq.  (16)  It  is  found  that 


fD 


(2.5) 


Comparison  with  the  data  of  Pig.  1  reveals  that  this  value  of  ffo  is 
similar  in  magnitude  to  that  actually  observed  during  test.  Test  data 
for  i  ■  1000  lie,  indicate  ffo  ■  2.1  (64/Njfc).  Prom  Eq.  (16)  for  the  mean 
velocity  it  ia  possible  to  define  an  effective  friction  factor. 

»e  .  u  (l  ♦  !!££)  (17) 


for  the  round  tube. 

Through  the  use  of  an  expression  for  viscosity  obtained  from 
the  kinetic  theory  it  is  possible  to  express  the  charge  rusher  in  terms 
of  the  Knuds en  number.  From  kinetic  theory  viscosity  can  be  expressed 


H  -  ax  N  m  c  A 


US) 
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where  ^  is  &  constant  equal  to  l/3  for  the  elementary  derivation  of 
viscocity.  Mobility  of  an  ion  can  be  expressed  as 


k.  5s_lA 

c 


(19) 


where  Sg  is  a  constant  depending  upon  the  theory  used.  Substituting 
for  and  K  in  the  charge  number,  one  obtains 

f2a.il  -  aa  JL  Sbl  (20) 

M  nS  n 


where  Ni/N  is  the  charge  density  ratio.  The  constant  *3  is  a  function 
of  ax  and  ae  and  has  a  value  of  the  order  of  unity.  It  can  be  seen 
that  the  Knudsen  number  acts  to  amplify  the  amount  of  charge  that  is 
present.  Since  Nrq  for  the  case  under  consideration  is  of  the  order  of 
10“°,  a  very  large  amplification  can  result.  The  foregoing  expression 
for  the  mean  velocity  is  analogous  to  an  expression  derived  by 
8teutzer.9 

An  extension  of  the  "induced  mean  body  force"  concept  is  very 
appealing  and  analyses  were  conducted  using  the  idea.10  Such  analyses 
provided  results  which  compared  well  with  the  test  data,  and  showed 
the  flattening  of  the  velocity  profiles  actually  observed.  If  the 
analyses  were  indeed  valid  it  raised  the  possibility  that  similar 
actions  could  be  observed  in  flat  channels  and  in  external  boundary 
layers  as  well. 

5.  Evaluation  of  Channel  Flow  Action 

Studies  were  then  initiated  in  two  separate  directions,  flat 
plate  channel  flow  and  external  boundary  layers  on  flat  plates.  The 
flat  channel  consisted  of  a  long  duct  5/8”  high  by  5"  wide  by  10  feet 
long.  Corona  wires  (0.004"  dia. )  were  stretched  lengthwise  and  were 
spaced  1/2"  apart.  Static  pressure  taps  were  located  along  the  length 
of  the  channel  to  measure  the  pressure  drop.  Figure  11  illustrates  the 
arrangement.  Data  obtained  on  the  pressure  drop  are  shown  in  Fig.  12 
where  it  can  be  seen  that  the  pressure  drop  is  linear  along  the  channel. 
Thus  it  does  not  appear  that  the  ion*  flow  action  in  the  channel  is  an 
end  effect.  Figure  13  shows  typical  friction  factor  vs.  Nr*  data  for 
the  flat  channel.  It  is  similar  to  that  found  for  the  round  pipe 
earlier. 


In  an  attempt  to  prove  or  disprove  the  hypothesized  existence 
of  an  induced  body  force ,  a  series  of  tests  were  run  to  determine  the 
current  distribution  along  the  length  of  the  flat  channel.  The  need 
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Fig.  12  -  Static  pressure  along  tbs  flat  duct 
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for  such  current  measurements  arises  from  a  consideration  of  the  elec¬ 
trical  field  equations.  The  current  involved  in  a  corona  discharge  le 
a  nonlinear  function  of  the  applied  voltage.  During  the  discharge 
phase  a  small  voltage  or  electric  field  change  can  lead  to  a  very  sig¬ 
nificant  change  in  the  corona  current.  Thus,  if  the  radial  or  trans¬ 
verse  E  field  is  varied  in  the  channel  locally  for  any  reason,  then  one 
would  anticipate  a  change  in  the  local  current.  Next  if  one  considers 
the  Maxwell  equation 


TOE  »  -  p 
dt 


for  the  case  with  B  small., 

TOE  -  0 


and  the  following  equation  is  significant. 

dEy  a  dEz 

Bz  By 

where  z  is  along  the  channel  and  y  is  transverse  as  shown  in  Fig.  -11. 

If  a  body  force  exists  in  the  channel  which  is  a  function  of  the  stressi 
velocity  along  the  channel,  that  is  Fz  »  f(w/k),  and  since  F*  *  pJSt, 
then  Ez  would  be  a  function  of  y,  and  BEy/Bz  could  be  non- zero.  Inis 
would  Imply  that  the  transverse  field,  Ey,  could  vary  along  the  channel. 

To  study  the  possible  variations  in  current,  the  flat  rhann^l 
was  modified  so  that  the  current  along  one  surface  could  be  monitored 
along  its  entire  length.  Figure  14  illustrates  the  modified  channel 
with  the  individual  conducting  silver  electrodes  painted  on  a  plastic 
sheet.  These  electrodes  can  be  seen  exposed  on  the  right  hand  side  of 
the  figure.  Two  traversing  pressure  probes  can  be  seen  on  the  left. 

In  conducting  the  tests  the  current  to  each  segment  was  measured  and 
the  local  current  density  determined.  The  results  are  plotted  in 
Fig.  15.  The  data  points  are  very  scattered.  If  a  significant  varia¬ 
tion  in  Ey  occurs  with  flow,  one  would  expect  to  see  a  distinct  change 
in  current  density  as  the  flow  velocity  is  increased  from  zero.  3tudy 
of  the  data  does  not  reveal  any  truly  marked  difference  in  current 
density.  The  data  do  indicate  however,  that  the  phenomena  are  not 
solely  "end  effects"  at  either  end  of  the  channel.  Thus  because  no 
clear-cut  change  in  current  density  was  found,  the  possiblity  of  the 
induced  body  force  was  seriously  called  into  question.  At  the  same 
time,  data  on  external  flow  testing,  which  will  be  described  subsequently. 
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15  -  Current  density  along  length  of  flat  channel 


_  .  _  • 


I 


also  indicated  that  a  significant  induced  streamwise  bod)'  force  prob¬ 
ably  did  not  occur. 

A  reconsideration  of  the  ion-flow  interaction  led  to  the 
consideration  that  the  changes  to  friction  factor  and  velocity  profiles 
in  the  channels  could  be  the  result  of  secondary  flows.  Accordingly 
a  two-pronged  effort  was  initiated,  involving  an  analytical  and  an 
experimental  study. 

The  analytical  approach  was  based  upon  flow  of  a  fluid 
between  concentric  cylinders  with  ions  moving  under  the  influence  of 
an  applied  field  in  a  radial  direction.  Because  of  the  collisions 
between  the  charges  and  the  fluid  an  electrical  body  force  in  a  radial 
direction  results.  Such  a  radial  force  may  be  considered  analogous  to 
the  centrifugal  force  in  Taylor  *s  problem  with  rotating  concentric 
cylinders  and  hence  it  was  hypothesized  that  a  pattern  of  internal 
vortices  could  possibly  be  induced.11* 12  The  coordinates  used  in  the 
analysis  are  shown  in  Fig.  l£.  The  analysis  started  with  the  Navier 
Stokes  and  included  the  radial  body  force  term  and  ion  diffusion  term. 
The  fluid  was  assumed  incompressible,  field  distortions  due  to  ion 
perturbations  small,  constant  conductivity  in  the  fluid,  fluid  proper¬ 
ties  unchanged  by  ionization,  and  that  the  perturbations  were  axlsym- 
metric.  After  a  perturbation  analysis,  the  applicable  equations  in  a 
nondime nsional  form  become 


+  o 

dr  r  Sz 


d  fSw_  Su"|  life- du"|  _  P  dg 

St  |_5r  Sz  J  L  r^J  |_dr  Sz  J  dz 


dn  .  Sir  ®i  n3  n 

+  u  ~  ”  -=■  V*n  m  0 

St  Sr  v 


Following  Taylor,  it  was  assumed  that 

u  ■  ux  (r)  cos  A  z  ejt 
w  ■  wx  (r)  sin  A  z  efft 
n  ■  nx  (r)  cos  A  z  eat 


. 

' 


where  X  and  o  are  the  noodimens ional  wave  number  end  amplification 
factor,  respectively.  With  substitution  of  these  assumed  forms  for 
the  perturbed  quantities  Eqs.  (24),  (25),  (26)  reduce  to 


'il  +  i-i-X* 

dr2  dr  r 


X2)  ux  -  }^Snx 


(30) 


ii*i±-X2 

dr2  r  dr 


') 


»i 


JL3£  ux 

Di  dr 


(31) 


The  boundary  conditions  of  the  Eqs.  (30)  and  (31)  are  that 
the  perturbation  quantities,  u,  v,  n,  induced  by  convection,  are  to 
vanish  at  each  of  the  bounding  cylinders.  Furthermore,  the  situation 
of  limiting  or  neutral  stability  is  given  by  <r  *  0.  The  problem  was 
next  restricted  to  a  limiting  case  where  the  distance  between  cylinders 
is  small  compared  to  inner  cylinder  radius,  and  an  analogy  with  Taylor 
problem  results.  In  this  case  the  terms  l/r  are  small  compared  to 
d/dr  and  the  above  equations  became 


(32) 


v  dn 
51  dr  Ul 


(33) 


The  noodlmenslonsl  independent  variable  r  was  changed  to  q  •  (r -r0)/l 
(r  dimensional)  and  Eqs.  (32)  and  (33)  were  combined,  for  the  case 
where  E  can  be  approximated  as  constant,  to  give 


(D2-a2)3  ux  ■  -  a2^! 


(34) 


where 


D  n  —  and  T  ■  -  —  X2£  £2 
dq  rj,3  dq 


is  the  electric  Taylor  nunfcer.  The  boundary  conditions  to  be  satisified 
by  ux  are 
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ui  -  DUX  •  (D2-4*)2Ui  ■  0  at  q  ■  ±  £ 


(35) 


1 
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0 
0 
0 
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0 
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0 
0 
a 
a 
a 
a 
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In  terms  of  dimensional  variables, 


T  takes  the  form 


14E  dfl~ 
^Dl  dr 


(36) 


It  can  be  shown  that  the  electric  Taylor  number,  T,  In  the  form  of 
Eq.  (36)  represents  the  ratio  of  the  d>  stabilizing  electric  force  to 
the  stabilizing  viscous  force.  The  solution  of  Eq.  (3*0  subject  to 
(35)  la  an  eigenvalue  problem  leading  to  a  relation  between  a  and 
T.  The  minimum  value  of  T  is  the  critical  Taylor  number  and  was  found 
to  be  1707.8  at  a  -  3.13.13 

Thus,  since  as  a  limiting  case  the  electrical  problem  reduced 
to  the  same  form  as  Taylor's  problem,  it  appeared  logical  to  expect 
some  type  of  internal  secondary  flows  to  form  within  the  channel. 

An  experimental  investigation  was  initiated  to  determine  If 
secondary  flown  were  present  in  the  channel.  Because  the  basic  flow 
velocities  were  of  the  order  of  a  few  feet  per  v  .  measurement  of 
the  very  small  secondary  flows  was  a  significsn*  problem.  Because  of 
the  intense  electrical  fields,  involved  use  of  hot  wire  anemoarntry  was 
ruled  out,  and  pitot  probes  were  insensitive  to  the  small  disturbances. 
Consequently,  recourse  was  made  to  a  technique  which  used  electrical 
probes  placed  in  the  fluid  stream.  The  arrangement  is  shewn  in  Fig.  17 
where  the  secondary  flow  as  hypothesized  is  shown.  The  two  probes 
shown  were  insulated  end  were  connect'.-  3  to  ground  potential  through  a 
resistor.  The  induced  secondary  flow  ^ould  assist  the  ion  motion  to 
the  end  surface  of  one  probe  and  oppose  ion  motion  to  the  second  probe. 
A  slight  difference  between  the  current  and  potential  of  the  two  probes 
should  result  if  some  type  of  vortex  pattern  exists  in  the  channel. 

This  potential  difference  was  applied  to  an  oscilloscope.  Typical 
results  are  shown  in  Figs.  8  and  19.  It  can  be  seen  that  fairly  regu¬ 
larly  spaced  oscillations  oo  occur  with  a  spacing  of  the  order  of  one 
diameter.  Thus,  it  appears  likely  that  same  secondary  motion  does 
exist  within  the  channel.  It  can  also  be  seen  that  as  the  Reynolds 
number  increases  the  regular  large  waves  disappear  as  indicated  in  the 
previous  friction  factor  tests,  Fig.  B. 

It  was  concluded  therefore,  that  the  ion-flow  interactions 
vritbin  the  cnannel  were  probably  the  result  of  the  creation  of  second¬ 
ary  flows  by  the  electric  wind  pressure  and  not  the  result  of  an 
induced  udal  body  force. 


n 


J 


31 


jlJL 


o  o  Q/ 

o  UtOO 


Pig.  17  -  Detection  of  induced  relodtie*  with 
elect  roetatic  probe* 
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Fig.  18  -  Probe  voltage  oscillations  Indicating  possible 
secondary  flow 
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6.  External  Flora 


Concurrent  with  the  internal  flow  study  of  the  nature  of  the 
ion-flow  interaction  in  the  channel!  a  program  was  initiated  to  etudy 
the  possible  existence  of  the  hypothesized  body  force  in  external 
boundary  layer  flow.  Analytical  and  experimental  investigations  were 
initiated.  The  analysis  was  based  upon  uaing  the  Prandtl  boundary 
layer  equation  for  a  flat  plate  including  a  stream  vise  body  force  of 
a  form  similar  to  Eqs.  (11)  and  (12). 11  The  equation  then  became 


v^+w^2-v^-£si 
by  di  dy2  pK 


(37) 


If  the  charge  density  is  assumed  constant  through  the  boundary  layer 
the  equation  takes  the  form 


r* 

by 


If* 


aw 


d*w 

SF 


(38) 


where 


This  equation  is  of  the  same  fora  as  Rossov's  equation  for  flow  of  a 
conducting  fluid  over  a  flat  plate  in  the  presenoe  of  a  transverse 
magnetic  field.14  Rossov's  solutions  indicate  that  the  boundary 
layer  is  distorted  by  the  magnetic  field  and  that  a  condition  of 
incipient  separation  can  be  reached. 

To  determine  the  possible  effects  experimentally,  a  heated 
flat  plate  was  mounted  in  a  flow  channel  and  the  effects  of  denotation 
and  a  transverse  field  were  studied. 19  The  plate  assembly  was  mounted 
in  a  Mach  Zehnder  interferometer.  Corona  vires  were  stretched  above 
the  surface  of  the  plate  in  a  atreasarise  direction  and  a  high  voltage 
applied  between  the  vires  and  the  plate.  The  plate  vas  heated  electri¬ 
cally  end  the  plate  surface  temperature  vas  monitored  by  the  use  of 
thermocouples  imbedded  in  the  plate  material.  The  physical  basis  for 
the  teats  vas  that  at  significant  stream  velocity  if  a  atreasarise  body 
fbroe  existed,  it  should  have  increasing  influence  on  the  boundary 
layer  as  the  ionization  vas  increased.  Since  the  thermal  boundary 
layer  is  a  good  measure  of  the  velocity  boundary  layer,  a  study  of  the 
fringe  patterns  obtained  should  provide  a  means  of  determining  the 
effects  of  the  ionization  on  the  external  boundary  layer.  The  experi¬ 
mental  arrangement  is  shown  in  Fig.  20. 
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The  most  significant  result  is  shown  in  Fig.  21  wherein  the 
ratio  of  the  Nusselt  number  with  ion  current  to  the  Ifueseit  number 
without  current  is  plotted  against  a  nondimensional  velocity.  This 
figure  shows  that  at  the  higher  stream  velocities  there  is  no  signifl- 
cant  influence  of  the  ion  currents.  At  low  velocities  the  effects  on 
heat  transfer  are  large,  but  this  can  be  expected  since  at  the  very 
low  velocities  the  electric  wind  actions  can  predominate  as  was  found 
from  the  previous  free  convection  tests  with  electrostatic  fields. 

These  tests  strongly  supported  the  conclusion  found  from  the 
studies  than  an  induced  axial  body  force  does  not  sema  to  exist. 

7.  External  Flows  -  Transition 

A  separate  effort  was  initiated  to  study  the  possible 
influence  of  an  electrostatic  field  (ions)  on  boundary  layer  transi¬ 
tion.-  A  flat  plate,  3*  x  6',  was  mounted  in  a  low  turbulence  tunnel 
and  corona  wires  were  suspended  above  the  plate  surface  as  shown  in 
Fig.  22.  Electric  fields  were  impressed  between  the  wires  and  the 
metal  surface  of  the  plate.  The  purpose  of  the  study  was  to  use  elec¬ 
trically  induced  perturbations  to  couple  with  the  natural  (Tbilasin- 
Schlichting)  oscillations  in  the  laminar  boundary  layer.  Corona  wire 
spacing  and  the  impressed  A-C  electric  field  component  were  selected 
to  match  the  predominant  wave  length  and  frequency  of  the  boundary 
layer  oscillations.  Boundary  layer  transition  was  located  by  means  of 
a  near  surface  impact  pressure  probe  that  was  traversed  slowly  from 
the  rear  of  the  plate  towards  the  leading  edge.  The  most  signiflcsnt 
results  are  shown  in  Fig.  23.  It  is  shown  that  transition  from  laminar 
to  turbulent  flow  as  indicated  by  the  low  point  in  the  curves  has  been 
shifted  backwards  along  the  plate  by  the  corona  discharge.  This  delayed 
transition  was  found  to  be  repeatable.  The  actual  mechanism  for  this 
transition  delay  is  still  the  subject  for  study.  It  is  interesting, 
but  no  effect  of  frequency  was  found  during  the  tests. 

8.  Flow  Attachment 

Because  the  energy  involved  in  electrostatic  actions  is 
usually  small,  it  is  necessary  to  consider  those  fluid  processes  where 
changes  to  the  flow  require  only  a  small  amount  of  energy  or  force. 

A  typical  situation  of  this  kind  occurs  in  the  fluid  attachment  as 
found  in  a  flip-flop  fluid  amplifier  element.  An  experimental  program 
was  initiated  to  determine  what  effects  electrostatic  actions  might 
have  on  attachment  and  the  flipping  action  of  a  flip-flop  unit  "  The 
experimental  equipment  for  a  large  exploratory  unit  is  shown  Fig.  2k, 
Air  passes  from  the  left-hand  side  through  a  plenum  chamber,  out  a 
nozzle  and  into  the  splitter  plate  region  to  the  right.  Two  side 
port  control  channels  are  shown  above  and  below  the  nozzle  exit  region. 
The  unit  was  five  inches  wide  with  a  nozzle  slot  height  of  0.1  inch. 
Several  electrostatic  arrangements  were  tried  but  the  most  interesting 
results  were  obtained  with  corona  wires  and  ground  electrodes  located 
in  each  control  channel.  In  operation,  a  voltage  waa  applied  between 
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the  corona  wire  and  ground  electrode  in  one  control  channel.  Hie 
electric  wind  then  iasued  from  the  aide  port  and  would  cause  the  main 
stream  jet  to  move  to  the  opposite  side.  When  the  corona  voltage  was 
switched  to  the  other  electrode  the  reverse  action  would  take  place. 

Thus  the  electric  wind  acted  similarly  to  the  side  port  air  input  for 
ordinary  flip-flop  fluid  elements.  The  power  required  to  flip  the 
main  stream  is  presented  in  Fig.  25.  Typical  downstream  velocity  pro¬ 
files  outward  from  each  aide  wall  cure  shown  in  Fig.  26. 

Because  of  the  success  in  achieving  flipping  action  with  the 
electric  wind,  a  brief  study  of  a  smaller  unit  was  accomplished. 18 
This  unit  had  a  main  nozzle  width  of  0.20  inch  and  a  height  of  l/k  inch. 
The  configuration  is  shown  in  Fig.  27.  The  study  included  several 
different  configurations  of  electrodes,  side  port  plenum  chamber  arrange¬ 
ments,  and  a  wide  range  of  flows.  In  general,  the  unit  operated  satis¬ 
factorily  with  typical  data  shown  in  Figs.  28  and  29.  Figure  28  illus¬ 
trates  that  powers  of  the  order  of  0.02  watt  are  necessary  to  flip  the 
mainstream  and  Fig.  29  indicates  a  flow  ratio  between  control  flow,  Qq, 
and  primary  flow,  Q,  of  the  order  of  0.125.  Attempts  were  made  to  deter¬ 
mine  the  maximum  frequency  at  which  the  unit  would  follow  the  input  elec¬ 
trical  pulaeo.  A  hot  wire  anemcneter  was  vised  to  determine  the  output 
flow  position,  but  unfortunately  considerable  difficulty  was  experienced 
with  the  unit  and  little  useful  data  were  obtained.  It  was  estimated 
that  the  upper  limit  of  the  device  based  upon  very  limited  testing  of 
the  order  10-20  cps. 


B.  CHANGE  OF  P^iASE 

The  possibility  of  using  electrostatic  fields  to  influence  the 
change  of  phase  of  a  fluid  is  particularly  intriguing.  During  change 
of  phase,  definite  surfaces  are  usually  present  where  surface  forces 
due  to  fields  can  be  induced.  There  is  also  the  large  difference  in 
the  number  of  molecules  per  unit  volume  between  the  phases  and  thus 
body  force  terms  could  also  be  active.  One  would  anticipate,  therefore, 
that  suitable  application  of  electrical  fields  could  alter  the  normal 
processes  during  the  change  of  phase  of  a  fluid.  Three  general  cases 
are  considered,  boiling  of  a  liquid,  condensation  of  a  vapor,  and  frost 
formation  upon  a  cold  surface. 

In  the  case  of  boiling,  very  interesting  work  has  been  done  by 
Bonjour  et  al. ,  Markels,  and  Choi. 19 »20* 21  Their  work  indicates  very 
large  increases  in  boiling  heat  transfer  rates  due  to  applied  electro¬ 
static  fields,  particularly  in  the  "burn-out"  region  of  the  boiling 
heat  transfer  curve.  Increases  in  heat  transfer  of  up  to  several 
hundred  precent  have  been  recorded  by  these  investigators.  They  have 
also  noted  a  complete  change  in  the  character  of  the  boiling  that  takes 
place. 
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1.  Condensation  of  Vapor 

The  condensation  of  vapors  in  the  presence  of  a  field  was 
studied  concurrently  and  independently  by  Choi  and  by  this  author.22* 

28  The  study  of  condensation  of  vapors  utilized  an  experimental 
arrangement  shown  in  Fig.  30.  Freon  113  was  boiled  by  means  of  the 
electrical  heaters  in  the  bottom  of  the  container.  The  cooled  copper 
plate  acting  as  the  condenser  surface  is  seen  on  the  side  of  the  con¬ 
tainer.  The  electrical  fields  were  achieved  by  placing  various  elec¬ 
trodes  directly  in  front  of  the  copper  plate.  Shown  in  the  sketch  is 
a  single  wire  electrode.  In  operation,  the  Freon  boiled,  and  then  the 
Freon  vapor  condensed  on  the  cooled  plate  and  ran  down  the  plate  back 
into  the  pool  of  Freon  at  the  bottom.  The  heat  balance  was  achieved 
by  very'  carefully  monitoring  the  mass  flow  rate  of  the  coolant  and  the 
change  in  its  temperature .  The  most  significant  results  were  achieved 
using  a  five-mesh  screen  grid  located  approximately  l/4  inch  ahead  of 
the  copper  plate.  The  results  are  shown  in  Fig.  31.  It  can  be  seen 
that  an  almost  threefold  increase  in  heat  transfer  occurs  with  the 
application  of  the  electrostatic  field.  The  line  drawn  represents  the 
Nusselt  theory  for  laminar  condensation.  The  power  required  to  accom¬ 
plish  the  increased  heat  transfer  is  indicated  in  Fig.  32  and  it  can 
be  seen  that  only  a  very  low  amount  of  energy  is  involved  in  obtaining 
the  increased  heat  transfer.  The  actual  condensation  process  with  the 
field  Imposed  is  very  interesting  to  observe.  The  condensation  test 
rig  is  brought  up  to  an  equilibrium  condition  and  the  voltage  is 
slowly  increased  between  the  screen  grid  and  the  cooled  copper  plate. 
Ho  observable  effect  of  the  field  is  noted  up  to  a  critical  value  of 
the  voltage.  Above  this  value  the  condensate  surface  first  becomes 
wavy  then  becomes  very  rough  in  appearance.  During  the  increased  con¬ 
densation  process,  approximately  as  much  condensate  is  observed  to  run 
off  the  screen  grid  an  runs  off  the  cooled  plate,  yet  no  Observable 
spray  is  seen  between  the  plate  and  the  grid.  Since  the  grid  remains 
at  vapor  temperature  it  is  obvious  the  condensation  process  does  not 
occur  on  the  grid. 

It  is  believed  that  the  mechanism  involved  is  similar  to  that 
of  colloid  electric  propulsion.24  Very  fine  charged  droplets  of  freon 
liquid  are  ejected  from  the  condensate  surface  and  are  accelerated  by 
the  electrical  field  across  the  gap  to  the  grid  electrode  where  they 
collect  and  run  off.  No  mist  is  seen  simply  because  the  droplets  are 
very  small. 

2.  Frost  Formation  and  Heat  Transfer  with  Fields 

A  simple  experimental  setup  was  used  to  investigate  the 
behavior  of  frost  formation  in  electric  fields.  The  setup  consisted 
of  three  parts:  the  test  plate,  made  of  a  l/2-inch  thick  sheet  of 
copper  and  having  a  frosting  surface  of  6  x  9  inches;  the  refrigerat¬ 
ing  system,  which  was  used  to  cool  the  test  plate  by  circulating  pre¬ 
cooled  refrigerant  through  the  back  of  the  plate;  and  the  high  voltage 


Fig.  30  -  Schematic  arrangement  of  experiments  of  condensation 
with  fields 
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:-lg.  32  -  Variation  of  condensation  heat  transfer  with  electrical 
power  input 


power  supply,  which  provided  the  strong  electric  fields.  The  test 
plate  was  installed  vertically  and  the  experiments  were  carried  out 
under  ambient  air  conditions. 

Preliminary  tests  were  conducted,  using  uniform  and  nonuniform 
fields  to  determine  the  possible  effects  of  electric  fields  on  frosting. 
The  uniform  field  was  achieved  by  applying  high  voltage  across  the  test 
plate  and  a  parallel  wire  screen.  A  nonuniform  field  was  established 
by  using  a  horizontal  corona  wire  parallel  to  the  centerline  of  the 
plate  with  the  wire  positive.  In  both  cases,  the  test  plate  was 
grounded.  The  results  of  the  preliminary  tests  can  be  summarized  as 
folio*  : 

1.  A  uniform  field  does  not  seem  to  have  any  observable 
effects  on  frosting  at  the  initial  stages  of  the 
process.  But  it  remains  to  be  determined  whether  a 
uniform  field  will  affect  the  growth  of  frost  at  a 
later  stage  when  loose  frost  begins  to  form. 

2.  A  nonuniform  field  characterized  by  a  corona  dis¬ 
charge  has  significant  effects  on  frost  formation. 

3.  For  the  case  of  wire-plane  electrode  configuration, 
the  frost  formed  a  strip  centered  along  the  line  of 
normal  projection  of  the  wire  on  the  plate.  In  the 
direction  transverse  to  the  strip,  the  frost  layer 
achieved  a  maximum  thickness  at  the  center  and 
thinned  out  in  both  directions  with  increasing  dis¬ 
tances  from  the  center.  Along  the  frost  strip  trans¬ 
verse  ridges  were  observed  to  appear.  These  may  be 
due  either  to  the  inherent  non- uniformity  of  the 
corona  or  to  an  electrohydrodynamic  instability 
induced  by  the  field.  Figures  33  and  34  show  effects 
on  frost. 

4.  The  e-fects  of  corona  discharges  on  frosting  can  be 
attributed  to  the  impingement  of  electric  wind 
induced  by  the  corona,  which  changes  the  flow  field 
from  free  convection  to  forced  convection. 

A  second  experimental  sequence  using  single  wire-plane  elec¬ 
trode  configuration  was  undertaken  to  study  the  flow  field  created  by 
the  impinging  electric  wind  and  to  provide  further  information  on 
electric  field  effects  on  heat  and  mass  transfer  processes  in  connec¬ 
tion  with  frost  formation.  The  mass  transfer  data  indicated  that  for 
a  test  period  of  five  minutes  an  increase  of  200$  in  the  total  frost 
accumulated  could  be  obtained  at  a  total  current  of  lOO^a,  as  compared 
to  the  case  of  no  field.  The  rate  of  increase  of  frost  deposition 
decreased  with  current,  suggesting  the  existence  of  an  upper  limit  on 
the  total  frost  versus  total  current  curves.  The  heat  transfer  data 
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also  showed  a  significant  increase  in  total  heat  transfer  rate  in  the 
presence  of  corona  discharge  with  an  increase  of  up  to  100^  at  150ua. 

The  velocity  distribution  of  the  electric  wind  was  measured 
by  means  of  a  total  head  pitot  traverse.  It  was  found  that  the  wind 
was  confined  to  a  narrow  region  within  one  inch  from  the  stagnation 
point  and  the  velocity  assumed  a  bell-shaped  distribution,  Fig.  35* 
Current  distribution  was  also  measured  underneath  the  corona  wire  by 
using  a  segmented  electrode  on  the  cooled  plate.  The  result  is  shown 
In  Fig.  36  for  various  values  of  applied  voltage.  Because  of  the 
similarity  of  the  two  curves,  velocity  and  current,  each  curve  was 
nondimens ionali zed  and  compared  assuming  a  centerline  value  of  unity. 
The  remarkable  agreement  is  shown  in  Fig.  37.  With  the  test  plate 
acting  as  a  cathode,  the  boundary  layer  velocity  profile  on  the  plate 
surface  was  measured  at  a  distance  two  inches  from  the  stagnation 
point  and  was  found  to  be  similar  to  that  of  a  wall  jet.  Because  of 
the  presence  of  the  strong  electric  field,  it  was  not  possible  to  mea¬ 
sure  the  velocity  distributions  in  the  stagnation  region  and  the  tran¬ 
sition  zone. 

The  flow  field  produced  by  an  impinging  electric  wind  and 
the  associated  heat  and  mass  transfer  processes  were  theoretically 
analyzed.  Comparisons  were  made  between  theoretical  predictions  and 
experimental  data.  Figures  38  and  39  show  that  the  correlation  was 
quite  good  for  the  case  of  frost  formation  with  applied  electrical 
fields  while  only  a  fair  degree  of  agreement  was  achieved  between  the 
heat  transfer  theory  and  data  with  fields  present.  In  suumary,  it 
appears  that  significant  changes  in  frost  formation  are  induced  by 
the  forced  convection  produced  by  the  electric  wind. 


C.  FLAME  HEAT  TRANSFER  WITH  FIELDS 

Another  possible  area  of  investigation  of  the  effects  of  electric 
fields  fall  in  the  area  of  the  effect  of  fields  on  flame  or  combustion 
heat  transfer.  During  combustion  of  relatively  high  degree  of  ioniza¬ 
tion  takes  place.  It  is  of  interest  to  determine  whether  this  nautral 
ionization  can  be  utilized  to  affect  the  heat  transfer  characteristics 
of  flames.  To  study  this,  flat  parallel  diffusion  flames  of  propane 
and  air  were  formed  in  a  flat  combustion  chamber;  the  wall  on  the  fuel 
side  served  as  an  anode  and  the  wall  on  the  air  side  was  the  cathode 
across  which  a  D.C.  voltage  of  3200  volts  was  impressed.25  Four  flames 
were  investigated;  three  flames  at  equal  velocities  with  various  pro¬ 
portions  of  propane  and  one  flame  at  a  higher  velocity.  The  leanest 
mixture  was  nonluminous  and  the  richest  mixture  very  luminous.  The 
combustion  chamber  had  segmented  anode  and  cathode  walls  and  was 
instrumented  to  measure  local  current  density,  local  heat  transfer 
rate,  electric  pressure,  exhaust  composition,  and  exhaust  temperature 
as  a  function  of  applied  voltage.  The  flames  could  be  viewed  through 
a  window  in  the  edge  wall  of  the  combustion  chamber. 
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Current  density  distribution  along  the  frosting  plate  surface 
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Figure  40  shows  the  experimental  setup  used.  The  upper  picture 
shows  the  parallel  sided  combustion  rig  with  the  many  heat  transfer 
coolant  loops  and  thermocouple  wires.  The  lower  picture  shows  the  gas 
analyzer  equipment  and  the  coolant  pump.  Figure  4l  shows  a  typical 
flame  under  various  levels  of  applied  electrical  field.  The  severe 
distortion  of  the  flame  is  clearly  evident. 

The  applied  electric  field  was  found  to  affect  the  flame  through 
three  different  mechanisms.  It  was  shown  that  a  gradient  in  current 
density  along  the  flame  was  responsible  for  an  induced  gas  motion 
away  from  the  flame  in  the  zone  containing  positive  ions.  Such  an 
induced  motion  affected  the  concentrations  of  fuel  and  air  in  the 
chamber  to  cause  the  flame  to  move  toward  the  cathode. 

Another  interaction  between  the  field  and  the  flame  was  the  onset 
of  a  flickering  condition  which,  with  increasing  voltage,  appeared  at 
upper  parts  of  the  flame  first.  A  parameter  for  describing  .’lame 
instability  was  set  up  which  indicated  that  conditions  near  the  satu¬ 
ration  voltage  of  the  flame  favored  instability.  A  stabilizing  para¬ 
meter  proportional  to  the  concentration  gradients  in  the  combustion 
chamber  was  defined  which  indicated  that  the  flame  gets  less  stable  in 
the  regions  further  downstream.  Onset  of  stability  was  detected  by  an 
unsteadiness  in  the  flame  current  observed  on  an  oscilloscope. 

A  third  mechanism  was  the  attachment  of  electrons  to  carbon 
particles  or  molecules  in  upper  parts  of  the  flame.  The  resulting 
electrlclally  Induced  drift  of  the  negative  ions  tended  to  reverse  the 
trend  of  induced  sotion  toward  the  cathode  at  lower  parts  of  the  flame. 

The  electric  interactions  with  the  flame  were  found  to  affect 
greatly  the  beat  transfer  rates,  especially  near  the  base  of  the  flsuse. 
At  the  highest  voltages  applied,  the  flame  was  observed  to  Impinge 
directly  on  the  cathode  with  increase  in  local  heat  transfer  over  the 
no  voltage  case  of  up  to  a  factor  of  6.  The  anode  heat  transfer  rate 
near  the  base  of  the  flame  was  found  to  be  reduced  up  to  85^.  At 
progressively  higher  parts  of  the  flame,  the  electric  effect  decreased. 
Typical  data  are  shown  in  Fig.  42.  It  is,  therefore,  concluded  that 
heat  transfer  rates  can  be  either  increased  or  decreased  in  a  controlled 
manner  by  the  application  of  fields. 

The  maxi  mum  current  density  from  the  flames  were  only  of  the  order 
of  either  microamperes  per  square  inch  of  flame  surface  and,  therefore, 
direct  effects  of  current  such  as  recombination  were  extremely  small 
cong)ared  to  the  heat  transfer  rates. 

Appreciable  decreases  in  the  amount  of  carbon  monoxide  and  unburnt 
hydrocarton  were  found  to  be  obtained  by  the  application  of  voltage; 
solid  carbon  formation  was  found  to  increase.  A  possible  explanation 
for  the  increase  in  solid  carbon  is  the  electron  attachment  to  carbon 
particles  in  the  pyiolysis  zone  with  a  subsequent  movement  away  from 
the  flame,  favoring  escape  of  the  particler  from  combustion. 
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There  were  indications  that  the  beat  transfer  was  convection  con* 
trolled,  but  a  possible  strong  mechanism  for  an  electric  effect  on 
luminous  flames  was  outlined.  It  nay  be  possible  to  control  the  emia- 
sivity  of  flames  through  the  influence  of  the  electric  field  on  the 
formation,  growth,  or  combustion  of  charged  particles  in  the  flame. 


D.  ION  FLOW  DIAGNOSTICS 

It  is  considered  possible  to  utilize  ions  injected  into  a  gaa 
stream  as  a  means  to  study  the  basic  nature  of  the  fluid  flow.  8ince 
in  the  case  of  air,  the  ions  would  be  either  0^  or  ions,  they  would 
tend  to  travel  with  the  fluid  readily  and  follow  the  particles  paths 
with  little  delay.  In  the  case  of  steady  flow  they  would  indicate  the 
actual  streamline  pattern.  As  an  illustration  of  their  use  in  this 
manner  one  has  only  to  look  at  the  study  of  the  secondary  flows  in  the 
channel  using  electrostatic  probes,  or  to  the  work  of  Donnelly  on 
secondary  flows.20  Jt  is  apparent  that  the  ion  currents  are  respon¬ 
sive  to  the  fluid  flow  situation  and  that  further  investigation  in 
this  direction  ia  warranted. 


1.  Velocity  Meter 

One  of  the  most  intriguing  applications  of  the  use  of  electro¬ 
static  effects  in  flow  diagnostics  is  the  air  velocity  meter  being 
studied  by  Professor  Durbin  of  Princeton. 27  In  this  device  the  air 
flow  to  be  measured  passes  through  a  cylinder.  At  the  axis  of  the 
cylinder  a  coaxial  electrode  is  placed  in  which  a  central  thin  disk  is 
placed  as  a  corona  discharge  electrode.  In  the  inner  wall  of  the 
cylinder  an  outer  electrode  is  located  which  has  a  finite  well  con¬ 
trolled  inteiral  resistance.  A  corona  discharge  is  initiated  between 
the  corona  electrode  end  the  cylinder  wall.  The  ions  created  tend  to 
move  radially  across  the  gap  towards  the  cylinder  wall.  The  air  flow 
moving  through  the  cylinder  carries  the  ions  downstream  and  they 
impinge  on  the  outer  wall  slightly  downstream  of  the  discharge  elec¬ 
trode.  When  this  happens,  the  potential  drop  along  the  wall  resistor 
is  now  different  for  the  forward  peurt  of  the  resistor  as  contrasted 
to  the  aft  portion.  By  means  of  a  suitable  external  bridge  and  an 
electrical  balancing  circuit,  the  device  will  indicate  the  velocity  of 
flow  which  caused  the  axial  drift  of  the  ions.  The  device  is  quite 
sensitive  and  it  is  believed  that  it  can  easily  measure  velocities 
down  to  a  few  feet  per  second.  An  eauviy  calibration  curve  is  shown 
in  Fig.  43. 

2.  External  Flow  Ion-Diagnostic  Studies 

In  an  attempt  to  extend  the  ion-diagnostic  techniques  to 
large  scale  flows,  a  series  of  tests  were  run  in  the  Army  Aeronautical 
Research  Laboratory's  7  x  10-Foot  Wind  Tunnel  at  Ames  Research  Center. 
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Initial  teats  were  conducted  on  a  six-inch  disaster  cylinder  three 
feet  in  length  counted  crossvise  to  the  sirstresa.  Ions  vere  created 
by  a  fine  vire  (0.0015"  dia.)  electrode  located  ahead  of  the  cylinder. 
Field  shaping  electrodes  were  placed  at  the  saw  forward  location  as 
the  corona  vire  but  above  anc*  below  the  vire.  Several  conducting 
electrode  strips  were  counted  on  the  cylinder  and  each  strip  was 
grounded  through  a  c canon  resistor  (12  cegoteu^  The  arrangement  is 
shown  in  Fig.  44.  The  electrical  output  of  ea  strip  was  monitored 
by  filtering  the  signal  developed  across  the  grounding  resistor  through 
a  Krohn-Hite  filter,  and  observing  the  wave  motion  on  an  oscilloscope. 
The  first  oscilloscope  picture,  Fig.  45,  shows  the  voltage  trace  at  a 
forward  position  on  the  cylinder  as  well  as  one  on  the  aft  part  of  the 
cylinder.  The  velocity  is  at  a  tunnel  q  of  one  (30*/Bec)*  The  regu¬ 
lar  oscillations  of  the  upper  traces  is  indicative  of  the  basic  flow 
about  the  cylinder.  These  oscillations  are  believed  to  be  caused  by 
potential  flow  oscillations  which  accompany  the  periodic  vortex  shed¬ 
ding  associated  with  the  Von  Karman  vortex  street.  The  lower  trace 
illustrates  the  flow  condition  in  the  aft  "dead  water"  region  of  the 
cylinder  and  it  is  very  clear  that  a  significant  difference  exists  in 
the  two  regions.  The  upper  traces  of  Fig.  46  i-lustrate  the  flow  at 
a  position  on  the  cylinder  just  forward  of  the  maximum  height  of  the 
cylinder,  and  the  lower  traces  show  the  flow  at  a  position  somewhat 
aft  of  the  maximum  position.  It  is  believed  that  the  upper  trace  still 
follows  the  overall  gross  oscillations  of  the  flow.  The  lower  trace 
showing  somewhat  erratic  motion  appear  to  indicate  a  separated  region 
of  flow  aft  of  the  maximum  height  of  the  cylinder.  Thus,  the  ion  tech¬ 
nique  seems  to  be  giving  proper  information  from  a  qualitative  stand¬ 
point  . 


Next  a  splitter  plate  was  mounted  directly  aft  of  the  cylinder 
at  the  centerline.  It  was  18  inches  long  and  fitted  with  an  aluminum 
plate  at  the  rear  to  provide  strength.  The  splitter  plate  was  used 
in  an  attempt  to  reduce  the  periodic  shedding  from  the  cylinder.  Ion- 
diagnostic  tests  were  run  with  this  configuration,  and  typical  results 
are  shown  in  Fig.  47.  The  upper  trace  represents  the  flow  near  the 
front  of  the  cylinder  and  the  lower  trace,  the  flow  in  the  aft  position 
near  the  juncture  of  the  splitter  plate  and  the  cylinder.  This  figure 
can  be  compared  with  Fig.  45  and  it  can  be  seen  that  the  oscillations 
axe  no  longer  as  cleanly  periodic  as  seen  previously  without  the 
splitter. 


Current  measurements  were  also  taken  at  the  aft  edge  of  the 
splitter  plate  in  an  attempt  to  determine  if  the  ion  measurement  tech¬ 
nique  could  be  effective  at  a  point  over  two  feet  away  from  the  corona 
wire.  In  Fig.  48  a  comparison  is  made  between  a  current  trace  from  an 
electrode  at  the  most  forward  side  of  the  cylinder,  and  a  trace  from 
the  rear  plate.  The  top  'Vide  trace"  is  from  the  front  of  the  cylinder 
and  the  lower  widely  oscillating  trace  is  from  the  aft  plate.  It  is 
clearly  evident  that  the  ion  technique  is  effective  at  significant 
distances  from  the  source.  It  should  be  noted  that  the  vertical  scale 
of  this  trace  is  reduced  by  one-half  from  previous  pictures. 
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Fig-  45  -  Oscilloscope  traces.  Current  q  ■  1-  Upper  trace  from 
forward  side  of  cylinder.  Lower  trace  frcm  rear  of  cylinder. 
Vertical  scale  2mv/cm 


Fig.  46  -  Oscilloscope  traces.  Current  3u&»  q  *  1.  Upper  trace  just  ahead 
of  maximum  cylinder  width  arid  lower  trace  just  aft.  Vertical 
scale  2mv/cm 


Next,  tests  were  run  with  a  flat  plate  mounted  in  the  7  x  10 
wind  tunnel  using  the  same  corona  wire  and  field  shaping  electrodes. 

The  arrangement  can  be  seen  in  Fig.  49.  Collecting  electrodes  were  pro¬ 
vided  on  the  surface  of  the  plate  by  using  conducting  silver  paint  as 
described  previously.  The  electrodes  were  strips  l/2  inch  wide,  l/3 
inch  apart  and  twelve  inches  long.  Each  was  grounded  through  a  12- 
megohm  resistor.  The  plate  was  aligned  parallel  to  the  stream  and  for 
a  tunnel  q  of  one  the  following  pictures  of  flow  oscillation  were 
obtained,  as  shown  in  Fig.  50.  The  top  trace  in  Fig.  50  illustrates 
the  oscillation  in  flow  four  inches  aft  of  the  leading  edge.  The  oscil¬ 
lations  may  indicate  the  possible  presence  of  a  separation  bubble  near 
the  leading  edge.  The  next  trace  from  a  strip  1-1/2  inches  aft  of  the 
first  trace,  although  it  shows  excursions,  does  not  show  significant 
oscillations.  The  remaining  traces  located  along  the  plate  to  a  posi¬ 
tion  approximately  30  inches  aft  of  the  leading  edge  are  relatively 
smooth  in  appearance. 

All  the  foregoing  traces  were  made  with  an  electrical  filter 
which  filtered  out  all  frequencies  above  15  cps.  Several  attempts  were 
made  to  locate  boundary  layer  characteristic  oscillations  of  the  Tollmein- 
Schlichting  type  or  the  turbulence  after  transition.  Unfortunately,  the 
high  frequency  noise  in  the  electrical  measuring  equipment  was  too  great 
and  no  useable  data  were  obtained  on  high  frequency  oscillations  in  the 
flow. 


Initial  tesus  were  also  conducted  on  an  NACA  0012  airfoil  in 
the  Army  Aeronautical  Research  laboratory  wind  tunnel.  A  corona  wire 
and  field  shaping  electrode  similar  to  that  used  previously  were  mounted 
in  front  of  the  airfoil.  Collecting  electrodes  were  located  on  the 
surface  of  the  blade  and  consisted  of  painted  silver  on  a  plastic  sheet 
in  a  manner  similar  to  that  used  on  the  flat  plate.  The  airfoil  used 
was  a  8-3/8  inch  chord  helicopter  tail  rotor  blade  and  was  mounted 
vertically  from  the  floor  of  the  wind  tunnel.  The  angle-of-attack  of 
the  airfoil  could  be  varied  and  the  corona  wire  installation  moved  with 
the  airfoil.  Initial  tests  at  a  tunnel  q  of  one,  and  an  angle  of  attack 
of  zero  degrees  indicated  erratic  electric  trace  behavior  for  traces 
near  the  nose  of  the  airfoil,  while  the  traces  along  the  upper  and  lower 
surface  were  generally  smooth.  The  erratic  trace  behavior  near  the  nose 
is  a  logical  subject  for  further  study.  When  the  airfoil  angle  of 
attack  was  raised  to  values  where  stall  should  occur,  the  traces  over 
the  "upper  surface"  leading  edge  showed  some  irregularity,  but  aft  of 
the  50^  chord  the  upper  surface  trace  oscillated  very  widely  and  errati¬ 
cally  indicated  the  presence  of  a  fully  separated  flow,  that  is,  ,,stall,, 
(Fig.  51). 


From  the  limited  work  accomplished  to  data  on  the  ion-flow 
diagnostic  technique,  it  appears  that  it  holds  great  promise  as  e. 
velocity  or  mass  flow  measuring  method,  it  holds  promise  in  indicating 
the  gross  nature  of  flow  disturbances  such  as  potential  flow  or  separa¬ 
tion  regions,  and  it  holds  promise  in  providing  the  integrated  history 
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Oscilloscope  traces  from  various  positions  along  the  surface 
of  the  flat  plate,  q  =  1,  i  =  5^a.  Upper  trace  is  four 
inches  back  of  leading  edges.  Remaining  trace  electrodes 
spaced  rearward  approximately  to  20  inches  aft  of  first 
trace  electrode 


: 

■ _ _  l 

of  flow  along  a  streamline.  With  instrumentation  available  to  date, 
however,  it  has  not  been  possible  to  study  the  details  of  the  boundary 
layer  flow  with  regard  to  frequency  oscillations  above  approximately 
25  cps. 


E.  POWER  GENERATION 

One  area  of  very  active  interest  in  electrofluidmechanics  is  that 
of  electrical  power  generation.  Although  this  investigator  has  not  per¬ 
sonally  conducted  studies  in  this  area,  a  brief  discussion  will  be  pre¬ 
sented  for  completeness .  Power  generation  with  electrostatic  is 
accomplished  in  a  manner  similar  to  the  familiar  Van  De  Graff  electro¬ 
static  high  voltage  generator.  Usually  a  closed  circuit  fluid  flow  is 
established  in  a  channel.  At  the  high  velocity  portion  of  the  channel 
ions  are  created  by  a  corona  discharge.  These  ions  are  carried  down¬ 
stream  by  the  high  velocity  flow  against  an  electrical  field.  The  flow 
thus  does  work  against  (upon)  the  charges  and  an  electrostatic  potential 
(work  per  unit  charge)  is  created.  Such  devices  usually  develop  very 
high  voltages,  low  currents,  and  are  D.C.  in  nature.  Their  primary 
advantage  over  the  MHD  power  generator*  lies  in  lower  temperature 
requirements  in  the  system.  The  group  under  von  Ohain  at  tile  Air  Force 
Aerospace  Research  Laboratories  has  been  very  active  in  this  area  of 
endeavor. 


IV.  SUMMARY 


The  entire  area  of  electro  fluidmechanics  offers  a  tremendous  chal¬ 
lenge  to  those  working  in  both  research  as  well  as  applied  fluid  mecha¬ 
nics.  The  range  of  possible  interactions  is  very  broad,  and  seems  to 
be  limited  only  by  the  limits  of  the  individual  investigator.  A  strin¬ 
gent  word  of  caution  is  suggested,  however,  to  those  who  may  wish  to 
pursue  research  in  this  field  of  endeavor.  The  energy  density  that  can 
be  created  by  an  electrostatic  field  is  very  small,  and  co-commitantly 
the  forces  that  can  be  developed  are  very  small.  If  onfe  attempts  to 
obtain  high  forces  with  a  field,  full  electrical  breakdown  will 
invariable  occur.  Secondly,  the  various  electrostatic  actions  tend  to 
occur  simultaneously.  For  example,  it  is  difficult  to  achieve  pure 
dielectrophoretic  action  without  some  pure  charge  phenomena  interacting 
at  the  same  time,  and  under  such  circumstances  it  is  only  too  easy  to 
misinterpret  the  results. 

It  is  believed,  however,  that  the  foregoing  summary  of  the  research 
that  has  been  conducted  does  indicate  that  real  opportunity  exists  for 
further  research  and  application  in  the  area  of  electrostatic-fluid 
interactions. 
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>  An  overall  review  of  studies  of  electrostatic  interactions  with  fluid  flows 
is  presented.  Review  of  prior  studies  involving  ions  with  free  convection,  and 
channel  flows  is  presented  along  with  prior  studies  of  field  effects  on  condensa¬ 
tion.  The  cause  of  large  changes  in  laminar  channel  flow  with  ionization  was 
studied  b-*th  analytically  and  experimentally.  The  changes  were  determined  to 
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the  action  of  an  impressed  E  field  were  found.  The  possible  use  of  ions  injected 
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a  very  logical  area  for  future  exploration,  t  - 
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and  flow  equations  is  attempted,  it  leads  to  a  complex  nonlinear  equa¬ 
tion.  Another  approach  based  upon  the  actual  experimental  problem  is 
possible.  It  can  be  assumed  that  the  low  velocity  flow  does  not  affect 
the  mechanism  of  the  corona  discharge  materially  and  consequently  the 
corona  charge  density  can  be  treated  independently  from  the  flow.  In 
addition,  because  of  the  length  of  the  channel  used  in  the  tests,  it 
can  also  be  assumed  that  the  transverse  charge  density  distribution  is 
essentially  constant  along  the  axis  of  the  channel.  Figure  10  illus¬ 
trates  the  coordinates  used. 


& 


The  current  equation  along  the  axis  of  the  channel  can  be 
written  as: 


Jz  =  pc  w  +  pc  K  Ez 


The  boundary  conditions  on  the  problem  require  that  the  current  path 
must  close  in  the  system  and  that  all  charge  that  leaves  the  wire  must 
reach  the  channel  wall.  This  condition  requires  that  no  net  current 
flows  out  of  the  channel,  iz  »  0  at  the  end  of  the  channel.  Hence 


J2  dA  a  0 


For  the  cylindrical  case  the  radial  current  density  is  given 


Jj*  a  ir/2  II  r  1 


Using  div  J  »  0,  for  cylindrical  coordinates. 


1  dr  Jr  .  d  Jz  n 
—  «  ■  ■■  4*  l  ■  —  U 

r  or  oz 


and  substituting  Jr  from  Eq.  (6) 


since  ir  must  be  a  constant. 
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Coordinate  system  for  round  tube 


Substituting  Eq.  (4)  into  Eq.  (5)  results  in 


w  +  pc  K  Ez)  dA 


0 


or 


i^pcEzdA=-i/Apc| 


dA 


(9) 

(10) 


It  can  be  recognized  that  each  integral  in  Eq.  (10)  represents  the  mean 
value  of  a  function.  Hence  we  can  write 


(Pc^z)m 


(Pcw)m 

K 


(11) 


and  noting  that  the  electrostatic  body  force  along  the  axis  of  the 
channel  is 


Fz  =»  PCEZ 


(12) 


Then  the  mean  value  of  pcEz  is  equal  to  the  mean  value  of  the  body 
force.  Substituting  into  the  fluid  flow  equation  for  channel  flow  with 
a  body  force  one  obtains 


p  vSf  -  |  (pcw)m  =  |£ 


(13) 


Non-dimension]  izing  this  equation  for  a  round  tube, 


d^  +  1  dw  (Pcv)m  N 


Pc  ■  NRe  S{ 


d£_ 


dtj  t]  drj 


Pcm  wm 


(14) 


where 


N 


Pc 


Pcm  K,  Pc*  -  . 


rj  -  r/R,  £ 


z/R,  W  .  jfc 


The  charge  number,  Npc,  is  a  new  dimensionless  parameter  which  char¬ 
acterizes  flow  with  ionization.  It  represents  the  ratio  of  the  charge 
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forces  to  the  viscous  force  in  the  flow.  Assuming  the  charge  density 
is  constant  one  can  solve  the  equation  to  get 


and  the  friction  factor  becomes 


(16) 


As  a  typical  case  consider  the  following  values,  K  =»  2  x  10**4 
m2/volt-sec,  R  *»  O.OI587  m,  i  a  10-3  amp,  channel  length,  i  a  1.7  meters, 
M  «  I.783  x  10"®  kg/m-sec,  then  pCm  =  2.j  x  109  charges/c. c.  The  charge 
number  becomes 


This  is  a  surprisingly  large  value.  Substituting  this  value  of  the 
charge  number  into  Eq.  (16)  it  is  found  that 


Comparison  with  the  data  of  Fig.  1  reveals  that  this  value  of  fp  is 
similar  in  magnitude  to  that  actually  observed  during  test.  Test  data 
for  i  *  IQ X)  jia,  indicate  ^  =  2,1  (64/NRe),  From  Eq,  (l6)  for  the  mean 
velocity  it  is  possible  to  define  an  effective  friction  factor. 

Me  =  M  (l  +  (17) 

for  the  round  tube. 

Through  the  use  of  an  expression  for  viscosity  obtained  from 
the  kinetic  theory  it  is  possible  to  express  the  charge  number  in  terms 
of  the  Knudsen  number.  From  kinetic  theory  viscosity  can  be  expressed 
as 

p  =  ax  N  m  c"  A  (18) 


CORRECTED 


20 


Fig.  46  -  Oscilloscope  traces.  Current  3n&,  q  =  1.  Upper  trace  just  ahead 
of  maximum  cylinder  width  arid  lower  trace  just  aft.  Vertical 
scale  2mv/cm 
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